On the global hydration kinetics of tricalcium silicate cement 
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We reconsider a number of measurements for the overall hydration kinetics of tricalcium 
silicate pastes having an initial water to cement weight ratio close to 0.5. We find that the time 
dependent ratio of hydrated and unhydrated silica mole numbers can be well characterized by 
two power-laws in time, x/(l — x) ~ (t/tx) ■ For early times t < t x we find an 'accelerated' 
hydration (%[> = 5/2) and for later times t > t x a 'deaccelerated' behavior (tjj = 1/2). The 
crossover time is estimated as t x ~ 16 hours. We interpret these results in terms of a global 
second order rate equation indicating that (a) hydrates catalyse the hydration process for 
t < t x , (b) they inhibit further hydration for t > t x and (c) the value of the associated 
second order rate constant is of mag nitude 6 • 10" 7 - 7 ■ 10" 6 liter mol' 1 s" 1 . We argue, by 
considering the hydration process actually being furnished as a diffusion limited precipitation 
that the exponents ip = 5/2 and tp = 1/2 directly indicate a preferentially 'plate' like hydrate 
microstructure. This is essentially in agreement with experimental observations of cellular 
hydrate microstructures for this class of materials. 

PACS number(s):81.30.Mh, 61.43.-j, 81.35.+k 



I. INTRODUCTION 



Heterogeneous solid-state transformations are of 
great importance and practical interest in industrial 
and technological applications. Perhaps most promi- 
nent examples arc transformations within alloys and 
steels playing a so crucial role for their mechanical and 
durabillity properties. Solid-state transformations are 
frequently accompanied by microstructural changes 
within the material in form of precipitating (segre- 
gating) phases from (solid) solutions. In its simplest 
form the transformation takes place between a major 
component acting as a solvent and an initially solvated 
phase. The precipitation process of the solvated phase 
is usually considered to happen in three main stages; 
phase nucleation, phase growth, and finally a coarsen- 
ing process leading to a homogenisation of microstruc- 
tural geometrical properties Further below we 
will employ basic concepts of solid-state transforma- 
tion in order to interpret the experimentally observed 
hydration kinetics of cement, See Section III. 



The massive introduction of hydraulic binders into 
daily engineering problems challenges a better under- 
standing on how these materials work. Beside this are 
cementing processes also of great importance for cer- 
tain geological problems, i.e., for the formation, prop- 
erties and behaviour of sedimentary bassins ]i| . Typi- 
cal hydraulic binders are plaster, cement, mortar, and 
concrete. In the following we will focus on the hy- 
dration of neat Portland cement (tricalcium silicate) 
more closely. 

Basic phenomenological aspects of cement hy- 
dration can be characterized as follows. Initially 
fine grained cement powder (here tricalcium silicate, 
C3S = Ca^SiO^, with grain diameters typically rang- 
ing between 5 jj,m and 50 p,m [Q) is mixed well with 
water. Tricalcium silicate is a well crystallized com- 
pound ( p = 3.21 g cm~ 3 ) and the employed pow- 



ders typically have specific surface areas of order of 
3 • 10 3 cm 2 g~ x ||. Rapidly after mixing the ce- 
ment particles start to dissolve. The principal re- 
action products are solvated ions (Ca 2+ , OH~ and 
H2SiO\~) diffusing within the solvent. The ion 
concentrations are bounded themselve by finite sol- 
ubility products above which hydrate phases start 
to precipitate from solution, preferable on surfaces 
of already existing hydrates. There are two associ- 
ated hydrates, a) the cement hydrate {C1.5SH2.5 = 
(CaO) U5 (Si0 2 ) (H 2 0) 2 .5) and b) the Portlandite 
(CH = Ca(OH)2), which compete for the common 
calcium and hydroxyl ions. The cement hydrate ( 
p = 2.35 g cm~ 3 ) is amorphous in contrast to the crys- 
talline Portlandite ( p = 2.24gcm -3 ). 

The process of cement dissolution, ion diffusion, and 
hydrate precipitation is usually referred as 'cement 
hydration'. The ion diffusion represents the phys- 
ical coupling between the chemical dissolution and 
precipitation reactions leading to a complex physico- 
chemical evolution of hydrate microstructure. 

In Fig. [I] we show a hydrated microstructure after 
40 hours of hydration time. Typical is the leaf/foil- 
like structure. 

In the following we shall propose on the basis of 
already published NMR measurements an empirical 
global reaction rate law for the hydration of cement 
pastes for initial water to cement weight ratios w/c = 
0.5. The observed values for the kinetics exponents 
can be understood - in parts - by classical solid-state 
transformation theory. 



II. EXPERIMENTAL RESULTS AND KINETIC 
RELATIONSHIPS 

The hydration of cement is spontaneous (exother- 
mic) and irreversible and can be characterized by the 
global net reaction HH, 



C 3 S (s) + AH 2 O w A C L5 Sff 2 , (s) + 1.5CiT (s) , (2.1) 

where k denotes an effective rate constant charac- 
terizing the solid to solid conversion rate. Because 
Eq. ( |2.l| ) describes a net reaction originating from 
vario us subprocesses the stoichiometric numbers in 
Eq. ( |2.l| ) are not related to the kinetic exponents of 
the overall reaction rate law. The kinetics of cement 
hydration has been experimentally studied to some 
extent by calorimetric and conductimetric measure- 
ments , b y X-ray diffraction J|,|| , by Raman spec- 
troscopy j[l()| , as well as by NMR spectroscopy |ll|-|l^| . 



A. NMR Measurements 

29 Si NMR spectroscopy has been used to exper- 
imentally determine the CSH growth rate ]Tl|-[T^|. 
The employed method (Single Pulse Excitation with 
Magic Angle Spinning (MAS/SPE)) is based on the 
fact that in the NMR spectrum the signal of the C3S 
silicon nuclei (monomer) is separated from those of 

The 



(dimer as well as trimer) mole fraction as Xb(P), i.e., 

XA{t) + X B {t) = 1. 



CSH silicon nuclei (dimers and trimers) |15 



distinction of dimer (silicon chain ends) and trimer 
(inside silicon chains) signals has been used to inves- 
tigate the polymerization process within the cement 
hydrate, see for example Ref. |lj. 

In order to enhance the signal/noise ratio one gen- 
erally needs to accumulate a great number of times 
the NMR signal. To achieve these accumulations and 
in order to obtain quantitative information about the 
populations of the different chemical species present in 
the specimen, one needs to strictly respect the relax- 
ation of each population between subsequent accumu- 
lations. We have selected from the literature results 
obtained taking this condition into account |Tl|-[l^|. 
NMR signal intensities, Ik, are in general proportional 
to the absolute number nk of excited silicon nuclei of 
species k within the sample, Ik = Tn^, with T being 
the constant of the measurement (species k denotes 
a silicon nuclei belonging to a monomer, dimer, or 
trimer silicon cluster). The proportionality constant 
is usually eliminated by considering relative signal in- 
tensities Ik/ Y2j Ij which correspond to mole fractions 
Xk = nk/ J2j n j ■ The constant Ylj n j defines the ab- 
solute silicon mass scale. The NMR measurements 
we refer to, provide plots of the time-dependent mole 
fractions for monomer, dimer and trimer silicon nuclei 
H-lia]. However, plots for mole fractions show, as 
will be demonstrated below, more complex algebraic 
behavior. 



B. Schematic Representation 

In order to get an idea about the global hydration 
kinetics of cement we consider the chemical net reac- 



tion (2T) in a more schematic representation, 



A^B, 



(2.2) 



where the species A and B refer to silica nuclei of 
C3S and CSH respectively. We denote the anhydrous 
C3S (monomer) mole fraction as XA{t) and the CSH 



The reader is reminded that Eq. (2.2) describes pre- 
cisely a solid-state transformation as mentioned in the 
introduction, however, there is one fundamental differ- 
ence in that the dimeric and trimeric silicon nuclei B 
are not initially solvated within the silicon monomers 
A. The B nuclei become rather formed within the 



liquid solution which does not appear in Eq. (2.2 



It is here the place to say something about the 
range of val idity, the similarities, and the differences 
of Eqs. (O) and @. 



Strictly speaking, the above net reactions can only 
hold, if the employed water to cement weight ratio 
is not too high, that is to say the number of sol- 
vated ions should be always orders of magnitudes 
lower than the number of molecules belonging to the 
solid state, i.e., Eqs. ( |2.1| ) and (2.2) are inconsistent 
in the limit of infinite dilution. A proximate condi- 
tion for this is a comparison of the cement to water 
concentration with the solubility of Portlandite, Sch, 
nc 3 s(to)/Vn 2 o(to) > S C h ~ 2 • 1(T 2 mol liter^ 1 , or 
w/c <C 200. The here considered experimental data 
certainly fulfill this condition, i.e., w/c ~ 0.5. 

Furthermore Eq. (|]l]) describes a sol id-liquid to 
solid transformation, however, Eq. ( |2.2| ) a solid to 
solid transformation. It is thus natural that both reac- 
tions describe very different things, if the initial water 
to cement ratio is so low that the water severely be- 
comes a limiting reactant for the hydration process 
(case of very thick pastes). To be more specific, the 
hydration following Eq. ( |2.l| ) cannot be completed for 
chemical reasons, if 4nc 3 s(io) > nH 2 o(to) or equiv- 
alently if w/c < (w/c)* = 0.31. It is, a priori, diffi- 
cult to judge whether the employed experimental wa- 
ter/cement ratio w/c = 0.5 is sufficiently large in or- 
der to exclude a systematic effect on hydration due to 
limiting water. However, as we will see further below, 
considering the kinetics of the 'deacceleration' period, 
water is very unlikely a chemically limiting reactant. 
It is, however, not possible to determine from the con- 
sidered measurements any water-specific kinetic expo- 
nent for the associated rate law of Eq. (2T). The de- 
termination of such a dependency, would require a set 
of experiments, conducted for different initial water 
to cement ratios. Similarly it is not possible to deter- 
mine kinetic exponents for Portlandite (CH) from the 
considered measurements, because ncsn is directly 
proportional to uch at every instant. In order to iso- 
late the kinetic influence of precipitated Portlandite 
on the overall hydration kinetics one needs to break 
up the direct proportionality between ucsh and uch ■ 
Preferable this could be done by examination of the 
hydration rate dependence for different initial CH ad- 
mixtures, i.e., ncnito) > 0. Bearing the foregoing 
caveats in mind, it should be obvious why we consider 
the schematic reaction (2.2). We are only in position 
to characterize the kinetic influence of reactants and 
products on the cement hydration. But already this 
limited information will be, in our opinion, informa- 
tive. 



C. Kinetic Relationships 



and 



In the following wc consider the ratio of silica mole 
fractions, 



f{t)=x B (t)/x A {t), 



(2.3) 



varying between (no hydrates at time t = 0) and 
oo (complete hydration at time t = oo). Figure || 
shows a double- logarithmic plot of f(t) versus the 
hydration time. This figure reveals several interest- 
ing facts about the hydration kinetics which will be 
addressed in this paper. We note that the consid- 
ered experimental data have been taken from three in- 
dependent publications obeying similar experimental 
conditions (room temperature, normal pressure, ge- 
ometry of specimens, and initial water/cement weight 
ratio w/c — 0.5), see jnHOl f° r details. The rela- 
tively narrow scatter of the data over order of magni- 
tudes demonstrates consistent and reproducible mea- 
surements. 

We observe, within in the scatter, two power laws 
separated by a sharply defined characteristic time 
i x = 16 hours with value f x = f(t x ) = 0.6. For 
the 'early' hydration period (t < t x ) we obtain, 



/(*) 



fx (tV, 



(2.4) 



with an exponent 7 = 2.5. Below 5 hours of hydra- 
tion time there are no quantitative NMR data avail- 
able, because of the bad signal/noise ratio for the 
CSH. Above t x = 16 hours the hydration dramat- 
ically slows down, 



/(*) = fx (f Y, 



(2.5) 



with S = 0.5. This second period takes place from 
16 hours to at least 8 • 10 3 hours of hydration time. 
Equations ( |2.4| ) and (2.1) define a continuous depen- 
dency in time. However, the first derivative of f(t) 
taken at t x is not continuous, which might appear 
unphysical. The discontinuity in df /dt arises, because 
we h ave assumed a zero crossover range in Eqs. (2.4) 
and (2.5) to obtain the best representation of experi- 
mental data. It is interesting to note, that an ansatz 
of the form 

t/t x = 2- 1 /™[(/// x ) Qm + (///x)' 3 " 1 ] 1 /™, (2.6) 

with proper chosen exponents a — 0.4 and j3 = 2.0 
gives a quite satisfactory representation of data for 
high integers of m, see Fig. |^ [To] . The best repre- 
sentation is found for m — > 00, which corresponds to 



Eqs. (gj) and (gj). 

In the following we propose relations between the 
above mentioned exponents 7 and 6 and the kinetic ex- 
ponents of an overall rate equation. The observed ex- 
ponents will allow for some general statements about 
the hydration mechanisms. 



We note, that from Eq. (2.3) the trivial relations 



x A (t) 



nA(t) 
n A (t Q ) 



! + /(*)' 



(2.7) 



x B (t) 



n B (t) 



f(t) 



n A (t ) ! + /(*)' 



(2.8) 



follow, where n A (to) characterizes the absolute silica 
mass scale. The overall rate equation for x A and x B 
has to be first order in time, 



1 



dt XA = dt XB 



(1 + f Y dt 



(2.9) 



One particular observation from Eqs. (2/7), (|2.8|) , and 
(2.9) is that the global hydration rate can be repre- 
sented as a product of powers of the mole fractions 
x A and x B , if 4zf follows a power-law in /. In such 
case the reaction is of overall order two and f(t) fol- 
lows a power- law in time, except for 4?f ~ / which 
yields an exponential dependency in time. We have 
already mentioned, considering the experimental re- 
sults in Fi g. 0, that fit) follows indeed two power- 
laws, Eqs. (|2.4"D and (^), seperated by a characteris- 
tic time t x . Thus it is possible to extract two rate laws 
from the measurements; one for 'early' times t < t x 
and another one for 'late' times t > t x . Consider 
fit) = f x ■ (t/t x ) where for times t < t x , rp = 7 and 
for times t > t x , ip = S, according to Eqs. (2.4) and 
( |2.5| ). The associated first order differential equation 

is (d/dt)f = VtxVx^/^ -15 ^- Inserting this into 
Eq. ( |2.9| ) and resorti ng w ith respect to factors x B and 
x A [Eqs. (2.8) and (|2.7D] we obtain the explicit rate 



equation in terms of mole fractions, 

" dt XA = dt XB = V ^>< /x X * 



(2.10) 



with tp = 2.5 for t < i x and ip = 0.5 for t > t x . 

Furthermore, by passing from mole fractions to con- 
centrations, it is possible to determine the approxi- 
mate effective rate constant (s) from the experimen- 
tally observed crossover point, the extracted expo- 
nents, and the initial water to cement weight ratio. It 
is convenient to consider first the initial concentration 
[A]\t — n A (to)/Vo of cement in the overall specimen 
volume V = Vc 3 s(to) + V H2 o(t ) @, 



1 



[A] I 



to 



= (1 + £C*£ . ^) WCaS pa 0.I8 liter mor 1 , 
Ph 2 o c 



with pc 3 s/PH 2 o 



3.21 and vc 3 s 



(2.11) 
7.2 ■ 



10 -2 liter mol~ x being the relative density and the 
molecular volume of cement respectively. The con- 
sidered water to cement weight ratio is w/c = 0.5. 
We find for the 'accelerated' period it < t x ), 



dt 



[A] 



dt 



[B] = k t<tx ■ [B]°- e [AY 



(2.12) 



with k t<tx = 2.5[A]\;H- 1 f° x A « 7 • 
10~ 6 liter mol~ x s _1 , and for the 'deaccelerated' pe- 
riod (t > t x ) 



with 



v t>t> 



[> [A 

Q.biAW^t- 1 /™ 



3.0 



(2.13) 
6 ■ 



III. INTERPRETATION 

Before we are going to present a possible expla- 
nation for the h ydra tion kinetics as manifested in 
Eqs. ( 2.12 ) and ( 2.13 ) we would like to give a brief 
account on another representation of the hydration 
kinetics being more widespread in cement literature, 
i.e., the degree of hydration a(t). The degree of hydra- 
tion is usually referred as the relative amount (mole 
fraction) of hydrated cemen t fig] , a(t) = £b(£) = 
1 — XA(t). Therefore Eq. (|2.8|) gives the relation- 
ship between a(t) and f(t) as defined in Eq. (|2.3| ). 
Typical experimental curves for the degree of hydra- 
tion exhibit sigmoidal shapes in linear representations. 
The observed inflection points, however, do in general 
not posses any particular significance for a change in 
chemical mechani sm. T his can be seen for example by 
assuming in Eq. fl2.10|) \& > 1 for all times. The de- 
gree of hydration will show an inflection point though 
there is only a single chemical mechanism (rate law) 
operative. Therefore it is not reliable to read off char- 
acteristic times from a(t) diagrams. While a(t) ap- 
proximates f{t) well if / <C 1 (early hydration) the 
discrepancy becomes very large for / 3> 1 (late hy- 
dration) . 

The foregoing remark has hopefully illustrated why 
we have avoided degree of hydration diagrams in our 
considerations. As another motivation for our ap- 
proach we present in Fig. || a so-called 'Avrami- 
Plot' for the hydrated silica amount, i.e., a test on 
the stretched exponential relationship xb = 1 — 
exp(—(t/r) k ). Such empirical relationships are fre- 
quently found for overall transformations Q|. It can 
be clearly seen from Fig. || that the hydration data 
cannot be described by such a r elatio n. 

The kinetic equations Eqs. (2.12) and ( 2.13| ) al- 
low to make somewhat more substantial statements 
about the global cement hydration mechanism(s). 
During the acceleration period already existing hy- 
drates catalyse the precipitation of new hydrat es (po s- 
itive kinetic exponent for the products in Eq. (2.12)). 
The growth of connected hydrate structures appears 
thus to be thcrmodynamically more favorable than an 
uncorrelated 'through solution' precipitation mecha- 
nism. 

On the other hand in the deacceleration period 
(t > t x ) the 'hydrate layers' surrounding the cement 
grains increasingly separate reacting anhydrous ce- 
ment and water and thus hinder/block further C3S 
dissolution. The hydrates are acting in the deaccel- 
erated period as inhib itors (negative exponent for the 
products in Eq. ( 2.13 ). The inversed role of hydration 
products during the accelerated and deaccelerated pe- 
riods appears to be experimentally evident from the 
foregoing considerations. 

A widespread assertion in cement literature is that 
'early' hydration is controlled by chemical kinetics 
whilst 'late' stage kinetics is being diffusion controlled. 
We agree with the later assumption that ion diffu- 
sion most probably represents the rate controlling 
step within the deacceleration period. This is in fact 
strongly indicated by the very low hydration rate at 
large times, see Fig. ||. 

However, there is no direct indication that the ac- 



celeration, kinetics is chemically limited. Obviously 
Eq. (2.12) cannot describe initial nucleation (approxi- 
mately within the first 15 min) because there exist no 
products at this times at all ([-B]|t = 0). On the other 
hand the early nucleation period is not accessible em- 
ploying the here considered experimental techniques, 
so there is no conflict in interpretation. One just has 
to keep in mind that all experimental data points as 
well as the thereof extracted power-laws are beyond 
the nucleation period. 

After the first few minutes of bringing cement and 
water into contact, the ions in solution rise their con- 
centrations far beyond the equilibrium solubilities, 
without precipitating at all. This can be understood 
by viewing the process of heterogeneous nucleation as 
overcoming a thermodynamical barrier (supersolubil- 
ity (7)). Being supersaturated the hydrate nucleation 
happens on the cement grain surfaces ]l9||l§| ]. At this 
instant a more or less significant part of the solution 
is strongly oversaturated with respect to the equilib- 
rium solubility. The further precipitation (growth) of 
the hydrates can thus be regarded to happen approx- 
imately within a spatial uniform oversaturated solu- 
tion (oversaturated interface layer). This is in so far 
of concern as the uniformity of 'initial conditions' is 
crucial to predict kinetic exponents for diffusion con- 
trolled reactions. It also supposes that the nucleation 
rate tends rapidly to zero as the hydrate microstruc- 
ture further develops. 

The hydrate microstructure has been experimen- 
tally classified for a water to cement ratio w/c = 0.47. 
Within the first 4 hours foil-like microstructure pre- 
cipitation is reported to happen radially away from 
the cement grains (typical dimensions < 0.5 \xm 
and x B < 10~ 2 ). Thereafter (< 24 hours) the 
formation of a gelatinous layer sourrouding the ce- 
ment grains have been observed (thickness ~ 0.5 /im 
and :eb(24 h) « 0.3). Also needle-like precipitates 
have been found. Finally after several days crum- 
pled interlocking foils are observed (comp. Fig. |l|). 
It has been also reported that the precipitated mi- 
crostructure morphologies are strongly influenced by 
the available interparticle spacings JlSj ] . There is 
presently no straightforward way to predict /calculate 
microstructural morphologies for such complex sys- 
tems, however, recently developed heterogeneous 
reaction-diffusion models have received considerable 
interest in this context [ pp[ . 

The growth of a precipitate is kinematically very 
complex. The precipitated geometry (microstructure) 
cannot be predescribed in general beyond the initial 
conditions but is rather a result of the associated inter- 
face dynamics. The most simplest case is governed by 
the growth of spherical precipitates of radius R from 
an initially supersaturated solution of concentration 
c. Let c = c(r) denote the particle concentration in 
the solvent, c s i the particle concentration within the 
precipitate at the interface, C£ S (R) the particle con- 
centration within the solvent at the interface and D 
the diffusion coefficient of particles in the solvent. One 
has from continuity of mass, 



dR 
~dt 



D 



dc 

cg s (R) \drJr=R 



(3.1) 



This is the equation of motion for the interface (in 
spherical co-ordinates). It establishes the direct pro- 
portionality between local growth rate and particle 
current density at the interface. In general C£ S de- 
pends, as a consequence of interfacial tension, on the 
(local) interface curvature. If the precipitate is not 
too small one often assumes the zero curvature limit, 
i.e., cjJR) w q s (oo). The most interesting point in 
Eq. ( |3.l[ ) is that the gradient of c contains 'global in- 
formation' about the interface because c(r) is the so- 
lution of the (stationary) diffusion equation, 



Ac(r) = 0, 



(3.2) 



= c and 



obtained under boundary conditions c(oo ) 
c(R) = ce a (oo). The solution of Eqs. (pTj) and (3.2) 
is the well known parabolic growth law 111], 



R*-Rl = 2D 5 -<*'(°°) ( t - 

Cat - Cfe(oo) 



(t-t ). 



(3.3) 



Similarly a parabolic solution is also obtained for the 
growth of a flat interface. More generally does the lo- 
cal growth rate dynamically depend on (at least) two 
competitive mechanisms (a) flattening of high curva- 
ture regions due to interfacial tension and (b) sharpen- 
ing of these regions due to preferential diffusive growth 
at these 'tips'. Numerical boundary integral methods 
have been developed in order to study the associated 
interface dynamics and instability, e. g. see Ref. [21]. 
For a pertubative treatment see for example Ref. [ 22 ] . 

Despite the complexity of involved microstructural 
transformations the natural question arises whether 
there exists at a given time a typical micro structure 
and a typical mode of growth within the system. 

Suppose the case of a vanishing nucleation rate dur- 
ing precipitate growth. If the precipitating structure 
grows geometrically in form of a plate (see the above 
experimental classification scheme) then the variation 
in mole fraction of precipitated phase in time is pre- 
dicted theoretically for xb -C 1 as ~ (t/r) 5 / 2 |jj 
with t being a characteristic timescale for the growth 
|23fl . The plate's rim grows at a constant rate while its 
thickness grows parabolically in time, explaining the 
exponent 5/2. For xb *C 1 one can replace the quan- 
tity f(t) by xs{t) in all foregoing considerations. One 
is lead to the conclusion that for early foil-lik e hy drate 
growt h 7 = 5/2 and (%p-l)/rp = 3/5 in Eqs. ( jgj ) and 
( 2.10| ) respectively. This is in agreement with the con- 
sidered measurements, see Fig. |^. 

What causes the observed slowing down of the hy- 
dration process? If a single foil would precipitate 
in a spatially infinite supersaturated solution there 
would be no obvious reason for a slowing down of 
the hydration process as there exists no characteristic 
length scale. However, the here considered case of ce- 
ment paste is an assembly of small anhydrous cement 
grains immersed in water. The mean free distance be- 
tween particles has to be considered as a typical length 
scale for the transport and for the precipitation pro- 
cess (xb <C 1). The size of the growing flakes can- 
not exceed this because of spatial hindrance. Hence 
there must exist a typical time scale t x at which the 
flakes change their mode of growth into thickening 
only. We interpret this typical time as the crossover 



) and (2.10) respectively, 
I 



time t x ~ 16 hours observed in Fig. Growth of 
flakes in the thickening only mode is expected to hap- 
pen parabolically in time for diffusion limited precip- 
itation reactions Therefore we predict the expo- 
nents of the deaccelerate d p eriod t o be S = 1/2 and 
(tp - l)/ip = -1 in Eqs. fT~ 
in agreement with the experimental data Fig. 

Apparently with the above assessments we have re- 
lated the experimentally observed kinetic exponents 
of the hydration products to microstructural informa- 
tion. Certainly there is no unique mapping between 
kinetics and geometry, bu t th is con stitu tes a complex 
question in terms of Eqs. (3.1) and (3.2) to be studied 
in future on its own right. 

So far we have restricted our considerations to the 
case of a water to cement weight ratio w/c = 0.5. The 
question arises whether the kinetic exponents are uni- 
versal and how the observed typical crossover quanti- 
ties t x and f x do depend on w/c. We show in Fig. f| 
experimental (replotted) data for three different w/c 
ratios. The data do not collapse, which is not so sur- 
prising because the typical interparticle spacing does 
depend on w/c. Qualitatively higher w/c ratios cor- 
respond to lower charateristic times t x and lower hy- 
drate 'amounts' f x . For early times we observe ki- 
netic exponents that do (apparently) depend on w/c. 
Interestingly are the kinetic exponents S for the deac- 
celeration period in all cases close to 1/2. However, we 
have not studied this in further detail because there 
are less experimental data available than for the stan- 
dard case w/c = 0.5, i.e., the hydration curves need 
to be experimentally reproducible for given w/c. 



IV. CONCLUSION 



We have reconsiderd NMR measurements on the 
overall hydration kinetics of tricalcium silicate pastes 
(w/c = 0.5) In Sec. II B we briefly discussed the con- 
ditions for w/c to be fulfilled in order to allow a mean- 
ingfull discussion in terms of a global net r eaction and 
its schematic counterpart. 



In Sec. II C we demon- 
strated that the time dependent ratio of hydrated and 
unhydrated silica mole numbers can be well character- 
ized by two power-laws in time, x/(l — x) ~ (i/i x ) ' ■ 
For early times t < t x we found an 'accelerated' hy- 
dration (tp = 5/2) and for later times t > t x a 'deac- 
celerated' behavior (tp — 1/2). The crossover time 
has been estimated as t x w 16 hours. We inter- 
preted these results in terms of a global second order 
rate equation indicating that (a) hydrates do catal- 
yse the hydration process for t < t x , (b) they do in- 
hibit hydration for t > t x and (c) the value of the 
associated second order rate constant is of magnitude 
6 • 10~ 7 — 7 • 10 -6 liter mol^ 1 s^ 1 . We have argued 
in Sec. Ill, by considering the hydration process actu- 
ally being furnished as a diffusion limited precipitation 
that the exponents -0 = 5/2 and -0 = 1/2 directly in- 
dicate a preferentially 'leaf like hydrate microstruc- 
ture. This argument was supported by experimen- 
tal observations of cellular hydrate microstructures for 
this class of materials. 
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FIG. 1. SEM micrograph of hydrated tricalcium silicate cement after 40 hours hydration time. The initial water to 
cement ratio was 0.5. Note the leaf-like hydrate microstructure (with courtesy of Institute Francais du Petrole). 
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FIG. 2. Double logarithmic plot of the measured mole ratio between hydrated and unhydrated silica f(t) = xb/ (1—xb) 
as a function of hydration time t employing NMR and TGA (Thermogravimetry). All measurements correspond to a 
water to cement weight ratio w/c — 0.5 conducte d u nder room temperature and normal pressure. The data are well 
represented by two power- laws: for t < t x by Eq. (BM (7 = 2.5) and for t > t x by Eq. (t^3) (5 = 0.5). The dashed line 
shows the corresponding numerical fit according to Eq. (|2.6|) employing t x = 16 hours, f(t x ) = 0.6, and m = 20. Data 
from (o) Ref. [n] (NMR), (x) Ref. (L2| (NMR), (□) Ref. [|J~(NMR) and (+) Ref. |TT| (TGA). 
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FIG. 3. 'Avrami'-Plot for the amount of hydrated silica (double logarithmic plot of the negative logarithm of unhy- 
drated silica amount — log(iBA) versus time t). If the hydration kinetics would follow a generalized Avarami-Johnson-Mehl 
law, i.e., xa = 1 — xb = exp( — (i/r) fc ), the data should lay on a single straight line of slope k. It can been seen that this 
is not the case. The data were taken from Ref. [12 (NMR) . 
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FIG. 4. Same representation as in Fig. ^ but for various water cement ratios employing XDR (x-ray diffraction) 
and Raman spectroscopy. The measurements were conducted at room temperature and normal pressure. The dashed 
line shows again the powerlaw fit corresponding to Fig. ^ for comparison purposes. Data from (□) Ref. [Hi] (Raman, 
w/c = 0.4), (o) Ref. H (XDR, w/c = 0.45) and (+) Ref. | (XDR, w/c = 0.7). Note that the data do not collapse, 
however, for large times they yield similar slopes ~ 1/2 compared to Fig. ^. 



